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. We analyze the physics potential of long baseline neutrino oscillation experiments 

\ planned for the coming ten years, where the main focus is the sensitivity limit to the 

' small mixing angle ^13. The discussed experiments include the conventional beam ex- 

. periments MINOS, ICARUS, and OPERA, which are under construction, the planned 

\ superbeam experiments J-PARC to Super-Kamiokande and NuMI off-axis, as well as 

Q^' new reactor experiments with near and far detectors, represented by the Double-Chooz 
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project. We perform a complete numerical simulation including systematics, correla- 



D , tions, and degeneracies on an equal footing for all experiments using the GLoBES 

' software. After discussing the improvement of our knowledge on the atmospheric pa- 

rameters ^23 and Aml]^ by these experiments, we investigate the potential to determine 
^ . 013 within the next ten years in detail. Furthermore, we show that under optimistic 

^ \ assumptions and for ^13 close to the current bound, even the next generation of ex- 

periments might provide some information on the Dirac CP phase and the type of the 
neutrino mass hierarchy. 



''Email: phuber@ph.tum.de 
''Email: lindner@ph.tum.de 
■^Email: rolinec@ph.tum.de 
■^Email: schwetz@ph.tum.de 
''Email: wwinter@ph.tum.de 



1 Introduction 



Within the last ten years a huge progress has been achieved in neutrino oscillation physics. 
In particular, the results of the atmospheric neutrino experiments [1-3] and the K2K accel- 
erator neutrino experiment [4] have demonstrated that atmospheric muon neutrinos oscillate 
predominately into tau neutrinos with a mixing angle close to maximal mixing. Further- 
more, solar neutrino experiments [5, 6] and the KamLAND reactor neutrino experiment [7] 
have established that the reduced flux of solar electron neutrinos is consistently understood 
by the so-called LMA-MSW solution [8]. Looking back at these exciting developments, it 
is tempting to extrapolate where we could stand in ten years from now with the experi- 
ments being under construction or planned. Certainly, neutrino physics will turn from the 
discovery era to the precision age, which however, will make this field by no means less 
exciting. The next major challenge will be the determination of the third, unknown mixing 
angle ^^13, which at present is only known to be small [9, 10]. Further important issues 
will be the determination of the neutrino mass hierarchy and, if ^13 turns out to be large 
enough, the Dirac CP phase. Three different classes of experiments are under discussion 
for the next generation of long-baseline oscillation experiments, which are able to address 
at least some of these topics: Conventional beam experiments, first-generation superbeams, 
and new reactor experiments with near and far detectors. In this study, we consider specific 
proposals for such experiments, which are under construction or in active preparation, and 
could deliver physics results within the next ten years. 

An already existing conventional beam experiment is the K2K experiment [4], which is 
sending a neutrino beam from the KEK accelerator to the Super-Kamiokande detector. 
This experiment has already confirmed the disappearance of as predicted by atmospheric 
neutrino data, and with more statistics it will slightly reduce the allowed range of the 
atmospheric mass splitting Am^i. In this study, we consider in detail the next generation of 
such conventional beam experiments, which are the MINOS experiment [11] in US, and the 
CERN to Gran Sasso (CNGS) experiments ICARUS [12] and OPERA [13]. These experiments 
are currently under construction and should easily obtain physics results within the next 
ten years, including five years of data taking. 

Moreover, we consider the subsequent generation of beam experiments, the so-called su- 
perbeam experiments. They use the same technology as conventional beams with sev- 
eral improvements. The most advanced superbeam proposals are the J-PARC to Super- 
Kamiokande experiment (JPARC-SK) [14] in Japan, and the NuMI off-axis experiment [15], 
using a neutrino beam produced at Fermilab in US. For these two experiments specific Let- 
ters of Intent exist and we use the setups discussed in there. JPARC-SK and NuMI could 
deliver important new results towards the end of the timescale considered in this work. 

Recently, there has been a lot of activity to investigate the potential of new reactor neutrino 
experiments [16]. It has been realized that the performance of previous experiments, such as 
CHOOZ [9, 10] or Palo Verde [17], can be significantly improved if a near detector is used to 
control systematics and if the statistics is increased [18-20]. A number of possible sites are 
discussed, including reactors in Brasil, China, France, Japan, Russia, Taiwan, and the US 
(see Ref. [16] for an extensive review). Among the discussed options are the KASKA project 
in Japan [19] at the Kashiwazaki-Kariwa power plant, several power plants in USA [21, 
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22] {e.g., Diablo Canyon in California or Braidwood in Illinois), and the Double-Chooz 
project [23] (D-Chooz), which is planned at the original CHOOZ site [10] in France. 

The particular selection of experiments considered in this study is determined by the re- 
quirement that results should be available within about ten years from now. This either 
requires that the experiments are already under construction (such as MINOS, ICARUS, and 
OPERA), or that specific proposals (Letters of Intent) including feasibility studies exist. 
From the current perspective, the only superbeam experiments fulfilling this requirement 
are the JPARC-SK and NuMI projects. Concerning reactors, we consider in this study the 
Double-Chooz project [23], since this proposal has the advantage that a lot of infrastructure 
from the first CHOOZ experiment can be re-used. In particular, the existence of the de- 
tector hall drastically reduces the required amount of civil engineering, which is considered 
to be time-critical for a future reactor experiment. Therefore, it seems rather likely that 
a medium size experiment can be built at the CHOOZ site within a few years and deliver 
physics results during the timescale considered here. We would like to stress that other 
reactor experiments of similar size, such as the KASKA project in Japan [19], would lead 
to results similar to Double-Chooz. To fully explore the potential of neutrino oscillation 
experiments at nuclear reactors, we furthermore consider an even larger reactor neutrino 
experiment (Reactor-ll). This could be especially interesting if a large value of 6*13 was found. 
Reactor-ll is the only exception for which we use an abstract setup, which could, in principle, 
be built at one of the sites mentioned above. For example, some projects discussed in the 
US, such as Diablo Canyon or Braidwood [22,24], are similar to our Reactor-ll setup. Such 
an experiment could be feasible within a timescale similar to the superbeam experiments, 
and could provide results at the end of the period considered in this work. Note that in this 
study, we do not consider oscillation experiments using a natural neutrino source, such as 
solar, atmospheric, or supernova neutrinos. 

The outline of the paper is as follows: After a brief description of the considered experiments 
in Section |21 we discuss the analysis methods and some analytical qualitative features of 
our results in Section El The main results of this study are given in Sections El El and [71 
First, in Sectional we investigate the improvement of the atmospheric parameters 623 and 
Amg^ from long-baseline experiments within ten years. Then we move to the discussion of 
the sin^ 2^13 sensitivity limit if no finite value of sin^ 2^13 can be established. We consider in 
SectionElthe conventional beam experiments MINOS, ICARUS, and OPERA. In SectionjHl 
we discuss the potential of reactor neutrino experiments to constrain sin^ 2^13, and we 
compare the final sin^20i3 bounds from the conventional beams, Double-Chooz, JPARC-SK, 
and NuMI. In Section [71 we investigate the assumption that sin^ 26*13 is large, and discuss 
what we could learn from the next generation of experiments on the Dirac CP phase and 
the type of the neutrino mass hierarchy. In this section, the Reactor-ll setup will become 
important. A summary of our results is given in Section |H1 In Appendix we describe 
in detail our simulation of MINOS, ICARUS, and OPERA. Furthermore, in Appendix iBl 
technical details of the reactor experiment analysis are given. Eventually, we present a 
thorough discussion of our definition of the sin^ 26*13 limit in Appendix O 
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2 Description of the considered experiments 

In this section, we discuss in detail the individual experiments considered in this work. The 
main characteristics of the used setups are summarized in Table E 

2.1 Conventional beam experiments 

Conventional beam experiments use an accelerator for neutrino production: A proton beam 
hits a target and produces a pion beam (with a contribution of kaons). The resulting pions 
mainly decay into muon neutrinos with some electron neutrino contamination. The far 
detector is usually located in the center of the beam. The primary goal of these beams is 
the improvement of the precision of the atmospheric oscillation parameters. In addition, an 
improvement of the CHOOZ limit for sin^ 2^13 is expected. For more details, see Ref. [11] 
for the MINOS experiment and Refs. [12,13] for the CNGS experiments. In addition, we 
describe our simulation in more detail in Appendix 1X1 

The neutrino beam for the MINOS experiment is produced at Fermilab. Protons with an 
energy of about 120 GeV hit a graphite target with an intended exposure of 3.7-10^° protons 
on target (pot) per year. A two-horn focusing system allows to direct the pions towards the 
Soudan mine where the magnetized iron far detector is located, which results in a baseline of 
735 km. The flavor content of the beam is, because of the decay characteristics of the pions, 
almost only z/^ with a contamination of approximately 1% z/g. The mean neutrino energy 
is at {Ey) ~ 3 GeV, which is small compared to the r-production threshold. The main 
purpose is to observe z/^ — > disappearance with high statistics, and thus to determine 
the "atmospheric" oscillation parameters. In addition, the z/^ — > Ve appearance channel will 
provide some information on sin^ 2^13. 

The CNGS beam is produced at CERN and directed towards the Gran Sasso Laboratory, 
where the ICARUS and OPERA detectors are located at a baseline of 732 km. The primary 
protons are accelerated in the SPS to 400 GeV, and the luminosity is planned to be 4.5 ■ 
10^^ pot y^^. Again the beam mainly contains z/^ with a small contamination of z/g at 
the level of 1%. The main difference to the NuMI beam is the higher neutrino energy. 
The mean energy is 17 GeV, well above the r-production threshold. Therefore, the CNGS 
experiments will be able to study the z/T--appearance in the z/^ — >■ v^- channel. Two far 
detectors with very different technologies designed for detection will be used for the 
CNGS experiment. The OPERA detector is an emulsion cloud chamber, whereas ICARUS 
is based on a liquid Argon TPC. In addition to the Vr detection, it is possible to identify 
electrons and muons in the OPERA and ICARUS detectors. This in addition allows to 
study the z/^ — > z/g appearance channel providing the main information on sin^2^i3, and 
the z/^ disappearance channel, which contributes signiflcantly to the determination of the 
atmospheric oscillation parameters. 

2.2 The first-generation superbeams JPARC-SK and NuMI 

Superbeams are based upon the technology of conventional beam experiments with some 
technical improvements. All superbeams use a near detector for a better control of the sys- 
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Label L 






oource 


Detector technology 


' ' ^Uet 


^run 


Conventional beam experiments: 










MINOS 735 km 


3GeV 


3.7- 


10^° pot/y 


Magn. iron calorim. 


5.4 kt 


5 yr 


ICARUS 732 km 


17GeV 


4.5- 


10^^ pot/y 


Liquid Argon TPC 


2.35 kt 


5yr 


OPERA 732 km 


17GeV 


4.5- 


10^9 pot/y 


Emul. cloud chamb. 


1.65 kt 


5yr 


Superb earns: 














JPARC-SK 295 km 


0.76 GeV 


1.0- 


10^1 pot/y 


Water Cherenkov 


22.5 kt 


5yr 


NuMI 812km 


2.22 GeV 


4.0- 


10^0 pot/y 


Low- Z- calorimeter 


50 kt 


5yr 


Reactor experiments: 












D-Chooz 1.05 km 


~ 4MeV 


2 X 4.25 GW 


Liquid Scintillator 


11.3t 


3yr 


Reactor- II 1.70 km 


~ 4MeV 




8GW 


Liquid Scintillator 


200 1 


5yr 



Table 1: The different classes of experiments and the considered setups. The table shows the label of 
the experiment, the baseline L, the mean neutrino energy {E^), the source power Psource (for beams: in 
protons on target per year, for reactors: in gigawatts of thermal reactor power), the detector technology, the 
fiducial detector mass mDct, and the running time trun- Note that most results are, to a first approximation, 
a function of the product of running time, detector mass, and source power. 

tematics and are aiming for higher target powers than the conventional beam experiments. 
In addition, the detectors are better optimized for the considered purpose. Since the pri- 
mary goal of supcrbcams is the sin^ 2^13 sensitivity, the z/^ z/g appearance channel is 
expected to provide the most interesting results. In order to reduce the irreducible fraction 
of Ve from the meson decays (which is also called "background") and the unwanted high- 
energy tail in the neutrino energy spectrum, one uses the Oj^- axis-technology [25] to produce 
a narrow-band beam, i.e., a neutrino beam with a sharply peaking energy spectrum. For 
this technology, the far detector is situated slightly off the beam axis. The simulation of 
the superbeams is performed as described in Ref . [26] ; here we give only a short summary. 

The J-PARC to Super-Kamiokande superbeam, which we further on call JPARC-SK,^ is 
supposed to have a target power of 0.77 MW with 10^"*^ pot per year [14]. It uses the Super- 
Kamiokande detector, a water Cherenkov detector with a fiducial mass of 22.5 kt at a 
basehne of L = 295 km and an off-axis angle of 2°. The Super-Kamiokande detector has 
excellent electron-muon separation and neutral current rejection capabilities. Since the mean 
neutrino energy is 0.76 GeV, quasi-elastic scattering is the dominant detection process. 

For the NuMI off- axis experiment [15], which we further on call NuMI, a low-Z-calorimeter 
with a fiducial mass of 50 kt is planned [28]. Because of the higher average neutrino energy 
of about 2.2 GeV, deep inelastic scattering is the dominant detection process. Thus, the 
hadronic fraction of the energy deposition is larger at these energies, which makes the low-Z- 
calorimeter the more efficient detector technology. For the baseline and off- axis angle, many 
configurations are under discussion. As it has been demonstrated in Rcfs. [26,29,30], a NuMI 
baseline significantly longer than 712 km increases the overall physics potential because of 
the larger contribution of matter effects. In this study, we use a baseline of 812 km and 

^The JPARC-SK setup considered in this work is the same as the setup labeled JHF-SK in previous 
publications [20,26,27]. 
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an off-axis angle of 0.72°, which corresponds to a location close to the proposed Ash River 
site, and to the longest possible baseline within the United States. The beam is supposed 
to have a target power of about 0.43 MW with 4.0 ■ 10^° pot per year. 

2.3 The reactor experiments Double-Chooz and Reactor-ll 

The key idea of the new proposed reactor experiments is the use of a near detector at a 
distance of few hundred meters away from the reactor core. If near and far detectors are 
built as identical as possible, systematic uncertainties related to the neutrino flux will cancel. 
In addition, detectors considerably larger than the CHOOZ detector are anticipated, which 
has, for example, been demonstrated to be feasible by KamLAND [7]. Except from these 
improvements, such a reactor experiment would be very similar to previous experiments, 
such as CHOOZ [10] or Palo Verde [17]. The basic principle is the detection of antineutrinos 
by the inverse /3-decay process, which are produced by /3-decay in a nuclear fission reactor. 
For details of our simulation of reactor neutrino experiments, see Ref. [20] and Appendix iBl 

For the Double-Chooz experiment, we assume a total number of 60 000 un-oscillated events 
in the far detector [23], which corresponds (for 100% detection efficiency) to the inte- 
grated luminosity of 288 1 ■ GW ■ yr, compared to the original CHOOZ experiment with 
12.25 1 ■ GW ■ yr leading to about 2 500 un-oscillated events [9]. The integrated luminosity 
is given as the product of thermal reactor power, running time, and detector mass. Note 
that, at least for a background- free measurement, one can scale the individual factors such 
that their product remains constant. The possibility to re-use the cavity of the original 
CHOOZ experiment is a striking feature of the Double-Chooz proposal, although it confines 
the far detector to a baseline of 1.05 km, which is slightly too short for the current best-fit 
value Am|i ~ 2 ■ lO'^eV^ 

If a positive signal for sin^ 26'i3 is found soon, i. e., sin^ 2^13 turns out to be large, it will be the 
primary objective to push the knowledge on sin^ 2^13 and 6cp with the next generations of 
experiments. From the initial measurements of superbeams, sin^ 2^13 and 6cp will be highly 
correlated (see Section In order to disentangle these parameters, some complementary 
information is needed. For this purpose, one can either use extensive antineutrino running 
at a beam experiment, or use an additional large reactor experiment to measure sin^ 26'i3 
precisely [20,31]. Because the antineutrino cross sections are much smaller than the neutrino 
cross sections, a superbeam experiment would have to run about three times longer in the 
antineutrino mode than in the neutrino mode in order to obtain comparable statistical 
information. Thus, a superbeam could not supply the necessary information within the 
anticipated timescale. We therefore suggest the large reactor experiment Reactor-ll from 
Ref. [20] at the optimal baseline of L = 1.7 km in order to demonstrate the combined 
potential of all experiments. It has 636 200 un-oscillated events, which corresponds to an 
integrated luminosity of 8 000 1 ■ GW ■ yr. Such a reactor experiment could, for example, be 
built at the Diablo Canyon or Braidwood power plants [22,24]. 
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3 Qualitative discussion and analysis methods 



In general, our calculations are done in the three flavor framework, where we use the standard 
parameterization U of the leptonic mixing matrix described by three mixing angles and one 
CP phase [32]. Our results are based on a full numerical simulation of the exact transition 
probabilities, and we also include Earth matter effects [8] because of the long baselines used 
for the NuMI beam. We take into account matter density uncertainties by imposing an 
error of 5% on the average matter density [33]. The probabilities are convoluted with the 
neutrino fluxes, detection cross sections, energy resolutions, and experimental efficiencies to 
calculate the event rates, which are the basis of the full statistical x^-analysis. We use all 
the information available, i.e., the appearance and disappearance channels, as well as the 
energy information. The simulation methods are described in the Appendices of Ref. [27]; 
for details of the conventional beam experiments, see also Appendix for the superbeam 
experiments Ref. [26], and for the the reactor experiments Ref. [20] and Appendix El All 
of the calculations are performed with the GLoBES software [34]. 

In order to obtain a qualitative analytical understanding of the effects, it is sufficient to use 
simplified expressions for the transition probabilities, which are obtained by expanding the 
probabilities in vacuum simultaneously in the mass hierarchy parameter a = Am'^i/Am'^i 
and the small mixing angle sin 2^13. The expression for the z/^ appearance probability 
up to second order in a and sin 26*13 is given by [35,36] 

~^ '^e) — sin^ 26113 sin^ 6*23 sin^ A 

=p a sin 2^13 sin Sep sin 2^12 sin 26^23 A sin^ A 

+ a sin 2^13 cos 6cp sin 2^12 sin 2^23 A cos A sin A 

+ cos^ ^23 sin^ 2^12 A^ (1) 

with A = Am|]^L/ {4:E,y). The sign of the second term is negative for neutrinos and positive 
for antineutrinos. The relative weight of each of the individual terms in Eq. is determined 
by the values of a and sin 26*13, which means that the superbeam performance is highly 
affected by the true values Am^i and Amg^ given by nature. Reactor experiments can be 
described by the corresponding expansion of the disappearance probability up to second 
order in sin 26*13 and a [19,20,36] 

1-Pee ^ sin^ 2^13 sin^ A + A^ cos^ ^13 sin^ 2^12. (2) 

The second term on the right-hand side of this equation is for sin^ 2^13 > 10~^ and close 
to the first atmospheric oscillation maximum relatively small compared to the first one, 
and can therefore be neglected in the relevant parameter space region. In principle, there 
are also terms of the order a sin^ 26*13 and higher orders in Eq. ^ . Though some of these 
terms could be of the order of the a^-term for large values of sin^ 26*13, they are, close to the 
atmospheric oscillation maximum, always suppressed compared to the sin^ 2^13-term by at 
least one order of a. Thus, the sin^ 26*i3-term carries the main information. 

From Eq. Q, it is obvious that a reactor experiment cannot access ^23, the mass hierarchy, 
or 6cp- In addition, the measurements of Aml-^^ would only be possible for large values 
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of sin^ 2^13 [20]. These parameters can be only measured by the z/^, z/^ — >• u^, and 

i^ti channels in beam experiments. However, comparing Eqs. and (j2)), one can 

easily see that reactor experiments should allow a "clean" and degenerate-free measurement 
of sin^ 26'i3 [19]. In contrast, the determination of sin^ 2^13 using the appearance channel 
in Eq. (jT)) is strongly affected by the more complicated parameter dependence of the oscil- 
lation probability, which leads to multi-parameter correlations [27] and to the (5, 6'i3) [37], 
sgn(Am3]^) [38], and (6*23, 7r/2 — 6*23) [39] degeneracies, i.e., an overall "eight-fold" degener- 
acy [40] . In the analysis, we take into account all of these degeneracies. Note however, that 
the (6'23, 7r/2 — 6*23) degeneracy is not present, since we always adopt for the true value of 6^23 
the current atmospheric best-fit value 6'23 = 7r/4. The proper treatment of correlations and 
degeneracies is of particular importance for the calculation of a sensitivity limit on sin^ 2^13. 
This issue is discussed in detail in Appendix m where we give also a precise definition of the 
sin^ 26'i3 sensitivity limit. In some cases we compare the actual sin^ 26'i3 sensitivity limit to 
the so-called (sin^ 2^i3)efT sensitivity limit, which includes only statistical and systematical 
errors (but no correlations and degeneracies). This limit corresponds roughly to the poten- 
tial of a given experiment to observe a positive signal, which is "parameterized" by some 
(unphysical) mixing parameter (sin^ 20i3)cff (see also Appendix [0 for a precise definition). 

If not otherwise stated, we use in the following for the "solar" and "atmospheric" parameters 
the current best-fit values with their 3 a-allowed ranges: 

I Am^, I = 2.Qtol ■ 10-3 eV^, sin^ 2^23 = I-0.15, 
Am^, = 7.0+?:^ ■ 10-^ eV^ sin^ 2^12 = 0.8l°f . (3) 

The numbers are taken from Refs. [41,42], which include the latest SNO salt solar neutrino 
data [6] and the results of the re-analysis of the Super-Kamiokande atmospheric neutrino 
data [2]. The interesting dependencies on the true parameter values are usually shown 
within the 3(T-allowed ranges. For the upper bound on sin^ 26'i3 at 90% CL (3cr) we use 

sin^ 2^13 < 0.14(0.25) , (4) 

obtained from the CHOOZ data [9] combined with global solar neutrino and KamLAND 
data at the best fit value Amg^ = 2 ■ 10"'^ eV^ [42]. In order to take into account relevant 
information from experiments not considered explicitly, we impose external input given by 
the la error on the respective parameters. This is mainly relevant for the "solar parameters" , 
where we assume that the ongoing KamLAND experiment will improve the errors down to 
a level of about 10% on each Am^i and sin 2^12 [43]. For the "atmospheric parameters" we 
assume as external input roughly the current error of 20% for |Am2j and 5% for sin2 2^23, 
which however, becomes irrelevant after about one year of data taking of the conventional 
beams, since then these parameters (especially |Am|^|) will be determined to a better 
precision from the experiments themselves. Furthermore, we assume a precision of 5% for 
|Am|^| for the separate analysis of the reactor experiments, since the conventional beams 
should supply results until then. However, it can be shown that the results would only 
marginally change for an error of 20% for lAmg^j. 

In general our results presented in the following depend on the assumed true values of the 
oscillation parameters. In particular they show a strong dependence on the true value of 
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Amgi, and therefore this dependence will be depicted in figures where appropriate. The 
^13 sensitivity limit obtained from Pg^ moreover also depends strongly on the true value of 
Am^i (see Figure IH below) . In principle also the variation of 6^12 plays a role. However, Pe^ 
depends only on the product of a ■ sin 26*12 up to second order in a as shown in Eq. ((H). 
Therefore a variation of the true value of 6*12 is equivalent to a rescaling of the true value 
of Am^i. The variation of the true value of 6*23 within the range given in Eq. (jS)) produces 
only slight changes in the results. In particular, those changes are much smaller than the 
ones caused by the variation of Aml^. Thus, in order to keep the presentation of our results 
concise, we do not explicitly discuss the dependence of the results on 623, and we always 
adopt the current best fit value 6*23 = 7r/4 for the true value. 

4 The measurements of Arri^^ and 623 

In this section, we investigate the ability of the conventional beam experiments and super- 
beams to measure the leading atmospheric parameters Am|^ and 6*23 . We do not include 
the reactor experiments in this discussion, since they are rather insensitive to Amli, and 
cannot access 6*23 at all. The measurement of these parameters is dominated by the i/^ — > z/^ 
disappearance channel in the beam experiments. 

In Figure we compare the predicted allowed regions for Amg^^ and sin^ 623 from the 
combined conventional beams (MINOS, ICARUS, OPERA), JPARC-SK, NuMI, and all beam 
experiments combined to the current allowed region from Super-Kamiokande atmospheric 
neutrino data. We show the fit-manifold section in the sin^ 6*23- Amg^^-plane (upper row), as 
well as the projection onto this plane (lower row). For a section, all oscillation parameters 
which are not shown are fixed at their true values, whereas for a projection the x^-function is 
minimized over these parameters. Therefore, the projection corresponds to the final result, 
since it includes the fact that the other fit parameters are not exactly known. In general, 
the x^-value becomes smaller by the minimization over the not shown fit parameters, which 
means that the allowed regions become larger. In Figure ^ the sgn(Am3]^)-degeneracy is 
not included, since it usually does not produce large effects in the disappearance channels. 
In addition, we use the true values sin^ 26*13 = 0.1 and 6cp = 0° in Figure [H Although 
the fit-manifold sections shown in the upper row of Figure Q depend to some extent on this 
choice, the effect for the final results of the disappearance channels is very small, i.e., the 
lower row of Figure ^ is hardly changed for sin^ 2^13 = 0. 

The first thing to learn from Figure ^ is that the precision on Amg^ will drastically improve 
during the next ten years, whereas our knowledge on 623 will be increased rather modestly. 
The combination of all the beam experiments will improve the current precision from the 
Super-Kamiokande atmospheric neutrino data [2] on sin^ 623 roughly by a factor of two, 
while the precision on Am|^ will be improved by an order of magnitude. Neither the three 
conventional beams combined nor NuMI will obtain a precision on 623 better than current 
Super-Kamiokande data, only JPARC-SK might improve the precision slightly. We note 
however, that the 623 accuracy of the long-baseline experiments strongly depends on the 
true value of Amg^, and it will be improved if Am'^i turns out to be larger than the current 
best- fit point. 
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Figure 1: The 90% CL (solid curves) and 3cr (dashed curves) allowed regions (2 d.o.f.) in the sin^ ^23- 
Am§i-plane for the combined conventional beams (MINOS, ICARUS, OPERA), JPARC-SK, NuMI, and aU 
beam experiments combined. For the true values of the oscillation parameters, we choose the current best-fit 
values from Eq. (j^J, a normal mass hierarchy, sin^ 2^13 = 0.1 and ^cp = 0. The upper row shows a section 
of the fit manifold (with the un-displayed oscillation parameters fixed at their true values), and the lower 
row shows the projection onto the sin^ 023-Am|]^-plane as the final result. The shaded regions correspond 
to the 90% CL allowed region from current atmospheric neutrino data [2]. 



In most cases, the correlations with the un-displayed oscillation parameters do not cause 
significant differences between the sections and projections in the upper and lower rows 
of Figure H Only for NuMI, the projection is affected by the multi-parameter correlation 
with sin^ 26*13 and 6cp- Since we do not assume additional knowledge about sin^ 26^13 for 
the individual experiments other than from their own appearance channels, the appearance 
channels can indirectly affect the Amg^^ or 623 measurement results. This can be understood 
in terms of the disappearance probability, which to leading order is given by [35,36] 

P,, = l-sin2 2^23sin2^^ + ... (5) 

where the dots refer to higher order terms in a = Am'^i/Am'^-^^ and 6*13, as well as products of 
these. Thus, the sin^ 26'13-precision, which comes from the appearance channels, is necessary 
to constrain the amplitude of the higher order terms in this equation which are proportional 
to 6*13. Since, however, 6*13 is strongly correlated with 6cp in the appearance channels, this 
two-parameter correlation can lead to multi-parameter correlations with 6*23 or Amg^^ in 
the disappearance channel through the higher order terms in Eq. (jSj). This explains the 
small differences between the section and projection plots in Figure ^ In addition, the 
measurement of sin^ 2^13 at NuMI is affected by matter effects, and hence, is somewhat 
different for the opposite sign of Am|^. Therefore, one can also expect a slightly different 
shape of the fit manifold for the sgn(Am3^)-degeneracy. Note that the initial asymmetry 
between sin^ 6*23 < 0.5 and sin^ 623 > 0.5 for NuMI is caused by its large matter effects. 
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Experiment / Combination 




^23 




sin 


2 a 


MINOS + OPERA + ICARUS 


O+0.34 
—0.18 


10" 


^eV^ 


(7r/4)^ 


-0.22 
-0.19 


0.5 


-un 91 
-0.18 


JPARC-SK 


9+0.15 
^-0.09 




c V 




0.5 


+0.13 
-0.10 


NuMI 


9+0.43 
^-0.07 


10- 


3eV2 


(71/4)^ 


-0.24 
-0.21 


0.5 


+0.23 
-0.20 


All beam experiments combined 


9+0.12 
^-0.06 


10- 


3eV2 


(vr/4)^ 


-0.13 
-0.10 


0.5 


+0.12 
-0.09 



Table 2: The expected allowed ranges (3(t, 1 d.o.f.) for the atmospheric oscillation parameters. For the 
true values of the oscillation parameters, we choose the current best-fit values, a normal mass hierarchy, 
sin^ 2^13 = 0.1, and Jcp = 0°. The impact of an inverted mass hierarchy, and different values for sin^ 20i3 
or (5cp on these final results is rather small. 

Eventually, one obtains the precision of the individual parameter Amg^^ or ^23 as projection of 
the lower row plots in FigureH(for one degree of freedom) onto the respective axis. In TableEl 
we show our prediction for the 3cT-allowed ranges of the atmospheric oscillation parameters 
from the conventional beam experiments and first generation superbeam experiments for 
one degree of freedom. 

5 Improved sin^ 2^13 bounds from conventional beams 

Let us now come to the crucial next step in neutrino oscillation physics: the determination 
of the small mixing angle ^13. We start this discussion by investigating the potential of 
the conventional beam experiments MINOS, ICARUS, and OPERA to improve the current 
bound on sin^ 26^13 . 

In Table El we show the signal and background event rates after one year of nominal opera- 
tion for each experiment (computed for sin^ 26'i3 = 0.1 and 5 = 0). Based on these numbers, 
one would expect that MINOS performs significantly better than ICARUS. However, TableEl 
only shows integrated event rates and does not include the energy dependence of signal ver- 
sus background event numbers. In the CNGS beam, the energy distribution of the intrinsic 
z/g-contamination is rather different from the energy distribution of the signal events. Thus, 
in a full analysis including energy information, the impact of the background is reduced. On 
the other hand, for the NuMI neutrino beam, the intrinsic i/g-contamination has an energy 
distribution which is much closer to the one of the signal events. Therefore, the impact of 
the background is relatively high. 

An important issue for the sin^ 26'i3 sensitivity limit from the conventional beams is the 
finally achieved integrated luminosity, which might differ significantly from the nominal 
value due to some unforeseen experimental circumstances. Therefore, we discuss the sin^ 26'i3 
sensitivity as a function of the integrated number of protons on target. In Figure El the 
sensitivity limits for MINOS, ICARUS, and OPERA are shown as a function of the luminosity. 
Note that since the CNGS experiments will be running simultaneously, we also show the 
combined ICARUS and OPERA sensitivity limit. In order to compare the achievable limits 
as a function of the running time, the dashed lines refer to the results after one, two, and 
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21.6 


12.2 
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0.33 


0.36 
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Table 3: The number of signal and background events after one year of nominal operation of MINOS, 
ICARUS, and OPERA. For the oscillation parameters, we use the current best-fit values with sin^ 26113 — 0.1, 
6cp — 0, and a normal mass hierarchy. 

five years of data taking with the nominal beam fluxes given in Refs. [12,44,45]. The 
lowest curves are obtained for the statistics limits only, whereas the highest curves are 
obtained after successively switching on systematics, correlations, and degeneracies. Thus, 
the actual sin^ 2^13 sensitivity limit in Figure |21 is given by the highest curves. The figure 
indicates that the CNGS experiments together can improve the CHOOZ bound after about 
one and a half years of running time, and MINOS after about two years. We note that 
the impact of systematics increases for MINOS with increasing luminosity, illustrating the 
typical background problem mentioned above. In Figure El we eventually summarize the 
sin^ 2^13 sensitivity after a total running time of five years for each experiment, assuming 
the true value of Amg^^ = 2 ■ 10^'^ eV^. One can directly read off this figure that the sin^ 2^13 
sensitivity limits of ICARUS and MINOS are very similar, and ICARUS and OPERA combined 
are slightly better than MINOS. 

Let us briefly compare our results to sin^ 26'i3 sensitivity limit calculations for MINOS, 
ICARUS, and OPERA existing in the literature. In the analysis of Ref. [46], the correlation 
with 5cp and the sign(Am||)-degeneracy are included, and hence these results should be 
compared with our final sensitivity limits, although we also include correlations with respect 
to all the other oscillation parameters. However, for the comparison, one has to take into 
account the different considered running times for MINOS (2 years vs. 5 years), as well as 
the difference in the chosen true value of |Am|^| (3.0 ■ 10"^ eV^ vs. 2.0 ■ 10"'^ eV^). In the 
analysis performed in Ref. [47], the correlation with (5cp and the sign(Am|]^)-degeneracy 
were not considered, while correlations with |Am|]^| were taken into account. Therefore, 
the results from that study should roughly be compared to our (sin^ 26'i3)efr limits. Again 
one has to take into account different assumptions about the running times and the true 
value for Am|]^ (5.0 • 10~^eV^ vs. 7.0 ■ 10~^eV^). Finally, in Appendix lA. 21 we demonstrate 
explicitly that our results are in excellent agreement with the ones obtained by the MINOS, 
ICARUS, and OPERA collaborations [12,44,45] if we use the same assumptions. 

A very interesting issue for the conventional beam experiments is the impact of the true 
value of Am2]^ on the sin^ 2^13 sensitivity. (The impact of the true value of Am|^ is dis- 
cussed in Sectional) One can easily see from Eq. (jT)) that the effect of 5cp increases with 
increasing a = Am^i/Amg^, which determines the amplitude of the second and third terms 
in this equation. Since the main contribution to the correlation part of the discussed figures 
comes from the correlation with 5cp (with some contribution of the uncertainty of the solar 
parameters), a larger Am^i causes a larger correlation bar. This can clearly by seen from 
Figure m which shows the combined potential of the conventional beams after five years of 
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Figure 2: The sin^ 29i3 sensitivity limit as function of the total number of protons on target at the 90% 
confidence level for MINOS, ICARUS, OPERA, and ICARUS and OPERA combined (5% flux uncertainty 
assumed). The dashed curves refer to the sensitivity limits after one, two, and five years of running. 
The lowest curves are obtained for the statistics limits only, whereas the highest curves are obtained after 
successively switching on systematics, correlations, and degeneracies, i.e., they correspond to the final 
sensitivity limits. The gray-shaded area shows the current sin^ 26i3 excluded region sin^ 2^13 > 0.14 at the 
90% CL [42] . For the true values of the oscillation parameters we use the current best-fit values Eq. Q and 
a normal mass hierarchy. 



running time (for each experiment) as a function of Amgi- In this figure the right edge of 
the blue band corresponds to the limit based only on statistical and systematical errors, i.e., 
the (sin^ 25i3)efr sensitivity limit. We find that the larger Am^^ is, the better becomes the 
systematics-based (sin^ 2^i3)efr sensitivity limit, and the worse becomes the final sensitivity 
limit on sin^ 26'i3. Since the LMA-II region is now disfavored by the latest solar neutrino 
and KamLAND data, Figure |3] demonstrates that the conventional beam experiments can 
definitively improve the current sin^ 26'13-bound. One may expect an improvement down to 
sin^ 26'i3 ^ 0.05 — 0.07 within the LMA-I allowed region, where the sin^ 26'i3 sensitivity limit 
at the current best-fit value is about sin^ 26'i3 < 0.06. 

Since the sin^ 2^13 sensitivity limit is expected to be sin^ 2^13 < 0.06 for MINOS, ICARUS, and 
OPERA combined (with five years running time for each experiment), a further improvement 
from the conventional beams seems to be unlikely. In addition, the systematics limitation, 
which can be clearly seen in Figure El demonstrates that a further increase of the luminosity 
would not lead to significantly better bounds on sin^26'i3. Therefore, one has to proceed to 



12 




sin^ 2^13 sensitivity limit 

Figure 3: The sin^ 26*13 sensitivity limit at the 90% confidence level after a running time of five years 
for the different experiments. The left edges of the bars are obtained for the statistics limits only, whereas 
the right edges are obtained after successively switching on systematics, correlations, and degeneracies, i.e., 
they correspond to the final sin^ 2^13 sensitivity limits. The gray-shaded area shows the current sin^ 2^13 
excluded region sin^ 26*13 0.14 at the 90% CL [42]. For the true values of the oscillation parameters, we 
use the current best-fit values Eq. @ and a normal mass hierarchy. 

the next generation of experiments to increase the sin^ 2613 sensitivity. Especially, the off- 
axis technology to suppress backgrounds and more optimized detectors could help to improve 
the performance. Amongst other experiments, we discuss the corresponding superbeams, 
which are using these improvements, in the next section. 

6 Further improvement of the sin 2^13 bound 

After the discussion of the conventional beams in the last section, we here discuss the final 
bound on sin^ 2^13 in ten years from now, if no finite value will be found (we will in the 
next section consider the case of a large 6*13). We first discuss in Section lUTTl the potential of 
a new reactor neutrino experiment, whereas we compare in Section IHUl the sin^ 2^13 limits 
from conventional beams, reactor experiments, and superbeams. 

6.1 Characteristics of reactor neutrino experiments 

In Figure El we show the sin^ 26'i3 sensitivity from reactor neutrino experiments as a func- 
tion of the integrated luminosity measured in t (fiducial far detector mass) x GW (thermal 
reactor power) x yr (time of data taking).^ We consider two options of the far detector base- 
line: Lfd = 1.05km, corresponding to the baseline of the CHOOZ site, and Lfd = 1.7km, 
which is optimized for values of Amg^ ~ (2 — 4) ■ 10~^eV^ [20]. A crucial parameter for 

^Note that we assume 100% detection efficiency in the far detector. For smaller efficiencies, one needs 
to re-scale the luminosity. 
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sin^ 2^13 sensitivity limit 



Figure 4: The sin^ 26*13 sensitivity limit at 90% CL for MINOS, ICARUS, and OPERA combined as 
function of the true value of Am^i (five years running time). The left curve is obtained for the statistics 
limit only, whereas the right curve is obtained after successively switching on systematics, correlations, and 
degeneracies, i.e., it corresponds to the final sin^ 2^13 sensitivity limit. The dark gray-shaded area shows 
the current sin^ 2^13 excluded region sin^ 2^13 > 0.14 at the 90% CL, and the hght gray-shaded area refers 
to the LMA-excluded region at 3ct, where the best-fit value is marked by the horizontal line [42]. For the 
true values of the un-displayed oscillation parameters we use the current best-fit values in Eq. ^ and a 
normal mass hierarchy. 

the sin^ 2^13 sensitivity is the uncertainty of the relative normahzation between the near 
and far detectors. We show the sensitivity for two representative values for this relative 
normalization uncertainty: First, Urei = 0.6% is a realistic value for two identical detec- 
tors [23]. Second, in order to illustrate the improvement of the performance of an reactor 
experiment with a reduced normalization error, we consider the very optimistic assumption 
of (Trei = 0.2%. Such a small value might be achievable with movable detectors, as discussed 
for some proposals in the US [24]. The shaded regions in Figure El correspond to the range of 
possible sensitivity limits for different assumptions of systematical errors and backgrounds. 
For the optimal case (lower curves), we only include the absolute flux and relative detector 
uncertainties (Tabs and (Jrei- For the worst limits (upper curves) we include in addition an 
error on the spectral shape o"shapc = 2%, the energy scale uncertainty (Xcai = 0.5%, and 
various backgrounds as discussed in Appendix iBl 

The first observation from Figure El is that the shaded regions in the left-hand panel are 
significantly wider than in the right-hand panel, which demonstrates that a reactor exper- 
iment at 1.05 km is more sensitive to systematical errors. This reflects the fact that the 
baseline of 1.7km is better optimized for the used value of [Amg^j = 2 ■ 10"'^ eV^, such 
that the oscillation minimum is well contained in the center of the observed energy range. 
In contrast, for the baseline of 1.05 km, the signal is shifted to the low energy edge of the 
spectrum. This implies that the interplay of background uncertainties, energy calibration, 
and shape error has a larger impact on the flnal sensitivity limit. 
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Figure 5: Luminosity scaling of the sin^ 20i3 sensitivity at the 90% CL. Here Am^^ = 2 • 10~^ eV^ is 
assumed to be known within 5%, Lnd = 0.15 km, and Lfd = 1.05 (1.7) km in the left (right) panel. The 
number of events in the near detector is fixed to 2.94 • 10^. We use Cabs = 2.5% and CTioI — 0.2% (0.6%) for 
the light (dark) shaded regions. The upper edge of each region is calculated for (Tshapc = 2%, deal = 0.5%, 
and backgrounds as given in Table El in Appendix ^ For the lower edges, we set (Tshapo = Ccai — and do 
not include backgrounds. The dots mark the Double-Chooz and Reactor-ll setups. 

However, from the left-hand panel, one finds that for experiments of the size of Double-Chooz, 
the impact of systematics is rather modest; the sin^ 2^13 sensitivity of 0.024 for normalization 
errors only deteriorates to 0.032 if all systematics errors and backgrounds are included. We 
conclude that the proposed Double-Chooz project is rather insensitive to systematical effects 
and will be able to provide a robust limit sin^ 29i3 ^ 0.032, although the far detector baseline 
is not optimized. In contrast, if one aims at higher luminosities, the systematics will have to 
be well under control at a non-optimal baseline such as at the CHOOZ site. In that case, it 
is saver to use a longer baseline. We note that the main limiting factor for large luminosities 
in the right-hand panel of Figure El is the error on a bin-to-bin uncorrelated background. 
Furthermore, comparing the light and dark shaded regions in that plot, it is obvious that 
a smaller relative normalization error will significantly improve the performance of a large 
experiment at 1.7 km, and will further reduce the impact of systematics and backgrounds. 
With the ambitious value of cXrei = 0.2%, sensitivity limits of sin^ 26'i3 < 7 ■ 10^^ could be 
obtained with a Reactor- ll-type experiment. 

Eventually, we have demonstrated that the Double-Chooz experiment could give a robust 
sin^2^i3 sensitivity limit. In fact, one can read off Figure El (D-Chooz-dot) that our as- 
sumptions about Double-Chooz are rather conservative. Since a Letter of Intent for this 
experiment is in preparation, we use it in the next subsection for a direct quantitative com- 
parison to the superbeams. However, as one can also learn from Figure luminosity and 
different systematics sources are important issues for a reactor experiment. Therefore, one 
should keep in mind that much better sin^ 26*13 sensitivity limits could be obtained from 
reactor experiments, such as the Reactor-ll setup. However, the exact final sensitivity limits 
will in these cases depend on many sources, which means that they are hardly predictable 
right now. 
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Figure 6: The sin^ 26*13 sensitivity limit at the 90% CL for MINOS, ICARUS, and OPERA combined, 
Double-Chooz, JPARC-SK, and NuMI. The left edges of the bars are obtained for the statistics limits only, 
whereas the right edges are obtained after successively switching on systematics, correlations, and degen- 
eracies, i.e., they correspond to the final sin^ 2^13 sensitivity limits. The gray-shaded region corresponds to 
the current sin^ 2^13 bound at 90% CL. For the true values of the oscillation parameters, we use the current 
best-fit values in Eq. ^ and a normal mass hierarchy. 

6.2 The sin^ 2^13 bound from different experiments in ten years from now 

Let us now assume that the conventional beam experiments MINOS, ICARUS, and OPERA 
have been running five years each, and that the Double-Chooz experiment has accumulated 
three years of data. In addition, we assume that the superbeam experiments JPARC-SK 
and NuMI have reached the integrated luminosities as given in Tabled (For earlier, more 
extensive discussions of the potential of superbeam experiments, we refer to Ref. [26].) 

In Figure El we show the sin^ 2^13 sensitivity for the considered experiments. The final 
sensitivity limit is obtained after successively switching on systematics, correlations, and 
degeneracies as the rightmost edge of the bars.^ Figure IHl demonstrates that the beam ex- 
periments are dominated by correlations and degeneracies, whereas the reactor experiments 
are dominated by systematics. It can be clearly seen that the (sin^ 2^i3)efj sensitivity limit 
(between systematics and correlation bar), or the precision of a combination of parameters 
leading to a positive signal, is much better for the superbeams than for the reactor experi- 
ments. Therefore, though the reactor experiments have a good potential to extract sin^ 26'i3 
directly, the superbeams results will in addition contain a lot of indirect information about 
6cp and the mass hierarchy, which might be resolved by the combination with complemen- 
tary information. We call this gain in the physics potential which goes beyond the simple 
addition of statistics for the combination of experiments "synergy" . In Section [3 we will 
discuss this further for the case if sin^ 2^13 turns out to be large. 

•^Note that earlier similar figures, such as in Refs. [20,26], are computed with different parameter values, 
which leads to changes of the final sensitivity limits. The largest of these changes come from the adjusted 
atmospheric best-fit values and NuMI parameters. 
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Another conclusion from Figure IHl is that it is very important to compare the sin^ 2^13 sen- 
sitivities of different experiments which are obtained with equal methods. In particular one 
clearly has to distinguish between the (sin^ 26'i3)efr sensitivity limit (between systematics and 
correlation bar) and the final sin^2^i3 sensitivity limit, including correlations and degenera- 
cies. For example, by accident the (sin^ 26'i3)efj sensitivity limit from the combined MINOS, 
ICARUS, and OPERA experiments is very close to the final sensitivity limit of JPARC-SK or 
NuMI. Thus, one may end up with two similar numbers, which however, refer to different 
quantities and are not comparable. 

A very important parameter for future sin^26'i3 measurements is the true value of Am|^, 
which currently is constrained to the interval 0.0011 eV^ < |Am|i| < 0.0032 eV^ at 3a [2]. 
From Figure!?! one can can easily see that the true value of Am|]^ strongly affects the sin^ 26'i3 
sensitivity limit. The left-hand plot in this figure demonstrates that for all experiments the 
sin^2^i3 sensitivity becomes worse for small values of lAmg^^l within the currently allowed 
range. However, since also the current sin^ 2^13 bound (dark-gray shaded region) is worse 
for small values of |Am3]^| than for large values, the relative improvement of the current 
sin^ 2^13 bound might be a more appropriate description of the experiment performance. 
This relative improvement as function of the true value of [Amg^l is shown in the right- 
hand plot of Figure!?! where a factor of unity corresponds to no improvement. From this 
plot, one can read off an improvement by a factor of two for the conventional beams, a factor 
of four for Double-Chooz, and a factor of six for the superbeams at the atmospheric best-fit 
value (vertical line). Nevertheless, the conventional beams might not improve the current 
bound at all for small values of lAmg^^l within the atmospheric allowed range, whereas any 
of the other experiments would improve the current bound at least by a factor of two. Thus, 
though the sin^ 2^13 sensitivity limit could be as be as large as sin^ 2^13 < 0.1 for small values 
of lAmgj^l for the superbeams or Double-Chooz, those experiments would still improve the 
current bound by a factor of two. 

7 Opportunities if sin^ 2^13 is just around the corner 

In Section !ni we have discussed how much the sin^ 2^13 bound could be improved if the 
true value of sin^ 26'i3 were zero. There are, however, very good theoretical reasons to 
expect sin^ 26*13 to be finite, such that the experiments under consideration could estab- 
lish sin^ 2^13 > 0. In this case, one could aim for the sin^ 26*13 precision, CP violation, 
CP precision measurements, and the mass hierarchy determination. Though it has been 
shown that CP and mass hierarchy measurements are very difficult for the first-generation 
superbeams and new reactor experiments [19,20,26,27,48], we will demonstrate in this sec- 
tion that we could still learn something about these parameters if sin^ 2^13 turns out to be 
large. In particular, we discuss the combination of the discussed experiments for the case 
sin^ 2^13 = 0.1. This would imply that a positive sin^ 2^13 signal could already be seen with 
the next generation of experiments. As discussed in Section !21 we assume here that a large 
reactor experiment Reactor-ll will be available at the end of the period under consideration 
to resolve the correlation between sin^ 26*13 and 5cp. We note again that similar results 
can be obtained by the superbeams in the antineutrino mode using higher target powers or 
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Figure 7: Left panel: The sin^ 2013 sensitivity limits at 90% CL from the experiments NuMI, JPARC-SK, 
Double-Chooz, and the combined conventional beams (MINOS, ICARUS, OPERA) as function of the true 
value of jATOgj^j. The dark-gray shaded region refers to the current sin^ 2^^i3 bound from CHOOZ and 
the solar experiments (90% CL) [42]. Right panel: The relative improvement of the sin^ 26*13 sensitivity 
limit with respect to the current bound from CHOOZ and solar experiments, where the dark-gray region 
corresponds to no improvement. The light-gray shaded regions in both panels refer to the atmospheric 
excluded regions (3cr), and the lines in the middle mark the current atmospheric best- fit value. 

detector upgrades.^ 

The superbeam appearance channels will lead to allowed regions in the sin^ 26'i3-5cp-plane, 
similar to the allowed regions for solar and atmospheric oscillation parameters from current 
data. We show the results of JPARC-SK, NuMI, and Reactor-ll for the true values sin^ 2^13 = 
0.1 and Sep = 90° in Figure |Hl and 6cp = —90° in Figure El For the right-most plots 
in these figures, we combine all experiments including the conventional beams MINOS, 
ICARUS, and OPERA, although they do not contribute significantly to the final result. 
Since we assume a normal mass hierarchy to generate the data, the best-fit is obtained 
by fitting with the normal hierarchy; the corresponding regions are shown by the black 
curves. The sgn(Am|^)-degenerate regions are obtained by fitting the data assuming an 
inverted hierarchy (gray curves). Thus, the best-fit and degenerate manifolds, which are 
disconnected in the six-dimensional parameter space, are shown in the same plots. Similar 
to Figure H we demonstrate the difference between a section of the fit manifold (upper rows) 
and a projection (lower rows) in these figures. 

As far as the measurement of sin^ 2^13 is concerned, any of the experiments in Figures |H1 
and El can establish sin^ 2^13 > for sin2 2^i3 = 0.1 at 3a. The sin^ 2^13 -precision can 

''in fact, one could already obtain some CP-conjugate information by running NuMI at L = 712km 
with antineutrinos only [48]. However, we do not consider an option with a very extensive a priori NuMI 
antineutrino running in this study, since the risk of this configuration is too high as long as sin^ 2^13 > is 
not established. 
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Figure 8: The 90% CL (solid curves) and 3a (dashed curves) allowed regions (2 d.o.f.) in the sin^ 26*13- 
(5cp-plane for the true values sin^ 26*13 = 0.1 and Sep = 90° for JPARC-SK, NuMI, Reactor-ll. The right-most 
plots are calculated for the shown experiments in combination with the conventional beams. For the true 
values of the un-displayed oscillation parameters, we choose the current best- fit values and a normal mass 
hierarchy. The black curves refer to the allowed regions for the normal mass hierarchy, whereas the gray 
curves refer to the sgn(Am§i)-degenerate solution (inverted hierarchy), where the projections of the minima 
onto the sin^ 20i3-5cp-plane are shown as diamonds (normal hierarchy) and triangles (inverted hierarchy). 
For the latter, the Ax^-value with respect to the best-fit point is also given. The upper row shows the fit 
manifold section (with the un-displayed oscillation parameters fixed at their true values) , and the lower row 
shows the projection onto the sin^ 20i3-(5cp-plane as the final result. 
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Figure 9: The same as Figure |H1 but for the true value (5cp 
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be read off from the figures as projection of the bands onto the sin^ 2^13-axis.^ The band 
structures of JPARC-SK and NuMI come from the CP phase dependency in Eq. (^. Because 
of the larger matter effects, the degenerate solution for NuMI is rather different from the 
best-fit solution, whereas it is very similar to the best-fit solution for JPARC-SK. For 
Reactor-ll, the sin^ 2^13-precision can be read off directly, since a reactor experiment is not 
affected by 6cp (see Eq. (j2I)), and the mass hierarchy has essentially no effect. Note that the 
treatment of the sgn(Am3^)-degeneracy in such a situation as shown in Figures |H1 and El is a 
matter of definition: One could either return two different intervals for normal and inverted 
mass hierarchies, or one could return the union of the two fit intervals as more condensed 
information. 

The figures show that for 5cp, none of the individual experiments can give any information, 
since no substantial fraction of antineutrino running is involved. However, there is some 
information on 6cp for the combination of all experiments, since the complementary infor- 
mation from the reactor experiment helps to resolve the superbeam bands. Note that the 
overall performance for the considered experiments (including the sgn(Am|^)-degeneracy) 
is usually better close to the true value Sqp = —90° than close to the true value Sep = 90°, 
since the degeneracy includes for 6cp = 90° very different values of 6cp far away from the 
best-fit manifold. This can, for example, be understood in terms of bi-rate graphs (c/., 
Refs. [38,48]). From a separate analysis of the CP precision, we find that one could exclude 
as much as up to 40% of all values of Sep at the 90% confidence level (1 d.o.f., close to 
Sep = —90°). However, if Sep turns out to be close to or vr, we find that one could not 
obtain any information on Sep- Furthermore, one can directly read off from Figures |H1 and IHl 
that CP violation measurements will not be possible with the considered experiments at 
the 90% confidence level (2 d.o.f. in these figures), since for the true values ^cp = ±90° 
corresponding to maximal CP violation, the projected allowed regions (even the best-fit 
solutions) include at least one of the CP conserving cases Sep G {0, 180°}. One can show 
that even for one degree of freedom, there is no CP violation sensitivity with the discussed 
experiments at the 90% confidence level. 

Another important issue for the next generation long-baseline experiments is the mass hi- 
erarchy determination. In Figures |H1 and IHl we give the Ax^-values for the minimum in the 
fit manifold corresponding to the sgn(Am|]^)-degenerate solution {i.e., for the inverted mass 
hierarchy) with respect to the best-fit minimum for the normal hierarchy (Ax^ = 0), which 
is the relevant number for the sensitivity to a normal mass hierarchy. Obviously, none of the 
individual experiments has a mass hierarchy sensitivity, but their combination has some. 
The mass hierarchy sensitivity becomes only possible because of the long NuMI baseline 
L = 812 km [26,29,30], since matter effects differ for the normal and inverted mass hierar- 
chies. Eventually, a NuMI baseline even longer than L = 812km could further improve the 
mass hierarchy sensitivity [26,48]. We note that the ability to identify the mass hierarchy 
strongly depends on the (unknown) true value of Sep- The mass hierarchy determination 
at the combined superbeams is close to the optimum for Sep = —90°, and close to the 

^Note that these figures are computed for two degrees of freedom, which means that the projections 
with one degree of freedom are slightly smaller. In fact, the 90% CL contour for two degrees of freedom 
(Ax^ = 4.61) is close to the 2a contour for one degree of freedom (Ax^ = 4.00). In particular, for the sake 
of comparison, we also use 2 d.o.f. for Reactor-ll, although it does not depend on 6cp- 
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minimum for Sqp = 90° [38,48]. In fact, one could have a better sensitivity to the normal 
mass hierarchy for 6cp = —90° (Ax^ = 4.9, see Figure EI) than for 6cp = 90° (Ax^ = 3.1, 
see Figure IHI) at the 90% confidence level (Ax^ = 2.71 for 1 d.o.f.). 

8 Summary and conclusions 

This study has focused on the future neutrino oscillation long-baseline experiments on a 
timescale of about ten years. The primary objective has been the search for sin^ 26'i3, but 
we have also analyzed the "atmospheric" parameters 623 and Aml^. The main selection 
criterion for the different experiments has been the availability of specific studies, such as 
LOIs or proposals, or that they are even being under construction. We assume that an 
experiment (including data taking and analysis) will only be feasible within the coming ten 
years, if it is already now actively being planned. The next long-baseline experiments will be 
the conventional beam experiments MINOS, ICARUS, and OPERA which are currently under 
construction. In addition, the JPARC-SK and NuMI superbeam experiments are under active 
consideration with existing proposals and will most likely provide results within the next 
ten years. Furthermore, new reactor neutrino experiments are actively being discussed. In 
this study, we have considered the Double-Chooz project, which will probably deliver results 
in the anticipated timescale, since infrastructure (such as the detector cavity) of the original 
CHOOZ experiment can be re-used. 

First, we have investigated the possible improvement of our knowledge on the leading atmo- 
spheric oscillation parameters. We have found that the conventional beams and superbeams 
will reduce the error on Amg^ by roughly an order of magnitude within the next ten years. 
The precision of 623 is dominated by JPARC-SK and will improve only by a factor of two 
(c/.. Table il). 

As the next important issue, we have investigated the potential of the conventional beams. 
I.e., MINOS, ICARUS, and OPERA, to improve the current sin^ 2^13 bound from CHOOZ 
and the solar experiments in a complete analysis taking into account correlations and de- 
generacies. Since the final luminosities of these experiments are not yet determined, we 
have discussed the results as function of the total number of protons on target. We have 
found that MINOS could improve the current bound after a running time of about two 
years, and ICARUS and OPERA combined after about one and a half years. In addition, 
we have discussed the maximal potential of all three conventional beams combined with 
a running time of five years each, leading to a final sensitivity limit of sin^ 2^13 < 0.061 
(all sensitivity limits at 90% CL). This final sensitivity limit includes correlations and de- 
generacies, which means that it reflects the experiment's ability to extract the parameter 
sin^ 2^13 from the appearance information. Since correlations and degeneracies could be 
reduced by later experiments, another interesting measure is the systematics-only sin^ 2^^i3 
limit for fixed oscillation parameters, i.e., the sensitivity limit to a specific combination of 
parameters, which we have called (sin^ 2^i3)eff. We have found a (sin^ 26'i3)efi sensitivity 
limit for the conventional beams of 0.026, illustrating that correlations and degeneracies 
have a rather large impact on the sin^ 2^13 limit from conventional beams. Note that it is 
important to compare different experiments with equal methods, which means that only 
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Current Beams D-Chooz JPARC-SK NuMI Reactor-ll Comb. 
sin^ 26*13 sensitivity limit (90% CL) 

sin2 26*13 0.14 0.061 0.032 0.023 0.024 (0.009) (0.009) 

(sin2 26*13 )eff 0.14 0.026 0.032 0.006 0.004 (0.009) (0.003) 
Allowed ranges for leading atmospheric parameters (3(t) 

|Am§i n+1.2 9+0.34 _ r,+0.15 9+0.43 _ 9+O.I2 

10-3 eV^ -0.9 ^-0.18 ^-0.09 ^-0.07 ^-0.06 

a ^7r\+0.20 /7rN+0.22 _ /7r\+0.13 itt\+0.2A _ /vrN+O.lS 

'^SS U^-0.20 14^-0.19 U^-0.10 V 4 .1-0.21 1 4 ^-0.10 

Measurements for large sin'^ 2^13 = 0.1 (90% CL) 

„;„2 9/1 Z. n 1+0.104 n -I +0.034 n 1 +0-067 p. +0.083 p. i +0.010 p, -,+0.010 

Sm ^t7l3 "-'.i_QQ52 U.1_QQ33 U.1_QQ34 U.1_QQ43 'J.l_Q_QQg U.i_QQQg 

5cp Combination can exclude up to 40% of all values 

CP violation No sensitivity to CP violation of any tested experiment or combination 
sgn(Am|i) Combination has sensitivity to normal mass hierarchy 



Table 4: Summary table of this study. The numbers which are printed boldface represent the best 
individual results within each row. For the true values of the oscillation parameters, we use the current 
best-fit values from Eq. (PJ and a normal mass hierarchy. The precisions for sin^ 26*13 do not include the 
sgn(Am|j)-degeneracy and are computed for the true value ^cp = 0. If one does not use Reactor-ll for the 
combination of all experiments, but Double-Chooz instead, one obtains the following values: 0.016 for the 
sin^ 2^13 limit, 0.003 for the (sin^ 29i^)cQ limit, and O.lj^Q Qji for the sin^ 20i3 precision. 



sin^ 26*13 or (sin^ 26*i3)efr limits should be compared with each other. In addition, it is in- 
teresting to observe that the final sin^ 2^13 sensitivity limit increases with increasing /S.m'^^ 
within the solar allowed region, whereas the (sin^ 26'i3)efr sensitivity limit decreases. This 
can be understood by the amplitude of the (5cp-terms which is proportional to Am\i. 

Furthermore, we have investigated the sin^ 2^13-limit obtainable by nuclear reactor exper- 
iments. A thorough analysis of the Double-Chooz configuration including systematics and 
backgrounds, demonstrates that a robust limit of sin^ 26*13 < 0.032 can be obtained in spite 
of the non-optimal baseline of 1.05 km. If one aims, however, to significantly higher lu- 
minosities than the 60 000 events anticipated by Double-Chooz, the systematics has to be 
well under control. In this case, a more optimized baseline of 1.7 km helps to reduce the 
impact of systematics and backgrounds, and limits of the order of sin^ 2^13 < 0.014 could 
be achievable. 

If in ten years from now no finite value is established, sin^ 2^13 bounds from the conventional 
beams (MINOS, ICARUS, OPERA), from reactor experiments, such as Double-Chooz, and 
from the superbeams JPARC-SK and NuMI will be available. We have demonstrated that the 
conventional beams could improve the current sin^ 2^13 bound by about a factor of two, the 
Double-Chooz experiment by about a factor of four, and the superbeams by about a factor 
of six. We have also shown that these results apply to a large range within the allowed 
interval for Am|i, since not only the experiment's potential decreases for small values of 
Am|i, but also the current sin^ 26*13 bound. For Am3i = 2 ■ 10~^eV^ we have found a 
final sin^ 26*13 sensitivity limit of sin^ 26*13 ^ 0.02 for the superbeams. Note that, though the 
Double-Chooz setup is not as good as the superbeams, its results are not affected by the true 
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value of Am2i within the solar-allowed range [20], which means that a reactor experiment 
is more robust with respect to the true parameter values. Moreover, because correlations 
and degeneracies do not effect the sin^ 29i3 limit from reactor experiments, the sin^ 29i3 and 
(sin^ 26'i3)efr limits are almost identical for Double-Chooz. In contrast, for the superbeams 
the (sin^ 2^i3)eff limit is nearly one order of magnitude smaller than the sin^26'i3 limit (c/.. 
Table EJ, demonstrating that correlations and degeneracies are crucial for them. 

In order to illustrate where we could stand in ten years from now if sin^ 2^13 were close to the 
current bound, we have also performed an analysis by assuming sin^ 26'i3 = 0.1. In this case, 
all the considered experiments will establish the finite value of sin^ 26*13 and measure it with 
a certain precision (c/.. Table Ej). In this situation, which is theoretically well motivated (c/.. 
Table 1 of Ref. [16]), one could even aim to learn something about 6cp and the neutrino 
mass hierarchy with the next generation of experiments. Since the results of superbeam 
experiments will lead to strong correlations between 6cp and sin^ 26*13, it is well known 
that complementary information is needed to disentangle these two parameters. One can 
either use extensive antineutrino running at the superbeams (which, however might not be 
possible at the time scale of ten years, because of the lower antineutrino cross sections), 
or a large reactor experiment to measure sin^ 2^13. In this study, we have demonstrated 
the potential for 6cp by assuming such a large reactor experiment at an ideal baseline of 
L = 1.7km, which we call Reactor-ll. Though such an experiment might not exactly fit 
into the discussed timescale, it might be realized soon thereafter. Possible sites for such an 
experiment are under investigation [16] (some proposals, which are, for example, discussed 
in the US, are similar to our Reactor-ll setup). Indeed, we find that in this optimal situation 
(sin^ 26*13 = 0.1), up to 40% of all possible values for 6cp could be excluded (90 % CL). 
This result, however, depends strongly on the true value of 6cp, and applies to maximal 
CP violation. For the case of CP conservation (true parameter value), however, nothing 
at all could be learned about 6cp- In either sensitivity to CP violation would 

not be achievable with the discussed experiments because of too low statistics. For the 
mass hierarchy determination, we have found that one would be sensitive to a normal mass 
hierarchy at the 90% confidence level, where the sensitivity to the inverted mass hierarchy 
would be somewhat worse [48]. Note that the sensitivity to the mass hierarchy is mainly 
determined by matter effects in NuMI, and could even be better for NuMI-baselines larger 
than 812 km. 

To summarize, from the current perspective neutrino oscillations will remain a very exciting 
field of research, and the experiments considered within the next ten years will significantly 
improve our knowledge. Eventually, these experiments could indeed restrict 6cp and deter- 
mine the neutrino mass hierarchy within the coming ten to fifteen years if sin^ 2^13 turns out 
to be sizeable. The remaining ambiguities could be resolved by the subsequent generation 
of experiments, such as superbeam upgrades, beta beams, or neutrino factories [29,49]. 
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A Simulation details of the conventional beam experiments 



In this appendix, we describe our simulations of the conventional beam experiments MINOS, 
ICARUS, and OPERA in greater detail. In Appendix lA.ll we give the numbers and refer- 
ences for the experimental parameters used, and in Appendix lA. 2! we demonstrate that our 
calculations reproduce the results of the simulations of the experimental collaborations to 
good accuracy. 

A.l Description of the experiments and experimental parameters 

The MINOS experiment will use both a near and far detector. The near detector allows to 
measure the neutrino flux and energy spectrum. In addition, other important characteristics, 
such as the initial Ve contamination of the un-oscillated neutrino beam can be extracted with 
good precision. Besides the smaller detector mass of Ikt, it is constructed as identical as 
possible to the far detector in order to suppress systematical uncertainties. The far detector 
is placed 713 m deep in a newly built cavern in the Soudan mine in order to suppress cosmic 
ray backgrounds. It is an octagonal, magnetized iron calorimeter with a diameter of 8 m, 
assembled of steel layers alternating with scintillator strips with an overall mass of 5.4 kt. 
The construction of the far detector was finished in spring of 2003 and it is now taking data 
on atmospheric neutrinos and muons. 

The mean energy of the neutrino beam produced at Fermilab can be varied between 3 and 
18 GeV. The beam is planned to start with the low energy configuration (PH21ow), with 
the peak neutrino energy at {Ej^) ~ 3 GeV. In our simulation we use the official PII21ow 
beam configuration [44], which means that we do not include a hadronic hose or different 
beam-plugs in the beam line setup. These modifications would lead to a better signal to 
background ratio. However, as discussed in Ref. [44], they affect the sensitivity limits to 
sin^ 2^13 only marginally. The NuMI PH21ow beam flux data, as well as the detection cross 
sections have been provided by Ref. [50]. We use 30 energy bins in the energy range between 
2 GeV and 6 GeV. In addition, the energy resolution is assumed to be aE = 0.15 ■ [11]. 
The NuMI beam will have a luminosity of 3.7 ■ lO^'^poty"^. In addition to the — > Ve 
appearance channel most relevant for the sin^ 26'i3 measurement, we include also the — > z/^ 
disappearance channel with an efficiency of 0.9, and we take into account that a fraction of 
0.05 of the neutral current background events will be misidentified as signal events. 

For the CNGS experiments, we use the fiux and cross sections from Ref. [51]. For both 
ICARUS and OPERA, we use an energy range between 1 and 30 GeV, which is divided into 
80 bins. For ICARUS, we assume an energy resolution of oe = 0.1 ■ Ej, [45], and for OPERA 
aE = 0.25 ■ Ey. The latter might be somewhat overestimated [13]. However, our limit 
at OPERA changes less than 5% for values of the energy resolution up to aE = 0.4 ■ E^. For 
the CNGS beam, we assume a nominal luminosity of 4.5 ■ 10^^ pot y^^. 

The original purpose of the CNGS experiments is the observation of — Vr appearance. 
The OPERA detector is an emulsion cloud chamber, and the extremely high granularity of 
the emulsion allows to detect the z/^ events directly by the so-called "kink", which comes 
from the semi-leptonic decay of the tauons. In order to reach a significant detector mass, 
the emulsion layers are separated by lead plates of 1 mm thickness. The total fiducial mass 
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of the detector will be 1.8 kt. However, during the extraction of the data, the detector 
mass will change as a function of time. Therefore, we use the time averaged fiducial mass 
of 1.65 kt for our analysis. A main challenge in the OPERA experiment is the automated 
scanning of the emulsions. The ICARUS detector uses a different approach: it is a liquid 
Argon TPC, which allows to reconstruct the three dimensional topology of an event with a 
spacial resolution of roughly 1 mm on an event by event basis. The detection is performed 
by a full kinematical analysis. The fiducial mass will be 2.35 kt. 

For the OPERA experiment, we include the information from the — > Vr channel by 
assuming an efficiency of 0.11, and a fraction of 3 ■ 10~^ of misidentified neutral current 
events. For the ICARUS experiment, we use an efficiency of 0.075 for this channel with a 
background fraction of 8.5-10"^ of the neutral current events [12]. Although the ICARUS and 
OPERA detectors are optimized to observe the decay properties of r-leptons, they also have 
very good abilities for muon identification, which allows to measure also disappearance. 
We therefore include the — > CC channel in both CNGS experiments, assuming a 
detection efficiency of 0.9 and taking into account a fraction of 0.05 of all neutral current 
events as background. As a matter of fact, the measurement of the atmospheric parameters 
also contributes to the sin^ 2^13 sensitivity limit, since correlation effects decrease with a 
higher precisions on the atmospheric parameters. Therefore, the sin^ 2^13 sensitivity at 
ICARUS and OPERA is considerably improved by including (besides the z/g channel) 

the — > Vr appearance and the z/^ disappearance channels in the fit. 

We have checked for all setups that the results do not depend significantly on the energy 
range, energy resolution, and bin size as long as the energy information is sufficient to 
distinguish the shape of the signal from the shape of the background. 

The sin^ 2^13 sensitivity of the different experiments is provided mainly by the information 
from the — > z/g appearance channel. Because of the small value of sin^26'i3, the number 
of Uf. CC events will be very small compared to the —>■ CC and NC events. 

Furthermore, the events from the intrinsic Ue component of the beam create a background 
to the oscillation signal. Thus, in our simulation, we consider as possible backgrounds: 
Beam z/g CC events, misidentified z/^ CC events, misidentified Ur CC events (mainly for 
CNGS), and misidentified NC events. We have calibrated the various background sources 
in our simulation carefully with respect to the information given in the literature. The 
corresponding references are for MINOS Table 3 of Ref. [44], for ICARUS Ref. [12], and for 
OPERA Table 4 of Ref. [45]. Using this information, we can reproduce with high accuracy 
the numbers of signal and background events provided by the experimental collaborations, 
which can be found in Table El 

A. 2 Reproduction of the analyses performed by the experimental collabora- 
tions 

In order to demonstrate the reliability and accuracy of our calculations, we use in this 
appendix the analysis techniques from Refs. [12,44,45], and compare our results with the 
ones in these references. For this purpose, we neglect all correlations and degeneracies, i.e., 
we set the solar mass splitting to zero, which also eliminates the solar mixing angle and CP 
effects. In addition, we fix the atmospheric mixing angle to n/A. Thus, the only remaining 
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Experiment Reference 


Signal 

I'e 


Background 
z/g ^ i/e l'^J,^ i>r NC Total 


MINOS NuMI-L-714 [44] 
ICARUS T600 proposal [12] 
OPERA Komatsu et al. [45] 


8.5 
51.0 
5.8 


5.6 3.9 3.0 27.2 39.7 
79.0 - 76.0 - 155.0 
18.0 1.0 4.6 5.2 28.8 



Table 5: The signal and background events for the three conventional beam experiments. The reference 
points are Am§i = 3.0 • 10-3eV^ sin^ 6*13 = 0.01 for MINOS, Am^i = 3.5 • lO^^gy^, sin^ 26*13 0.058 
for ICARUS, Ato§i = 2.5 • 10-3eV^ sin^ 26I13 = 0.058 for OPERA, and sin^ 26I23 = 1, Am|i = sin^ 26li2 = 
5qp = in all three cases. The nominal exposures are lOkty (MINOS), 20kty (ICARUS), and 8.25 kty 
(OPERA). Note that these numbers are different from the ones in TableOl since different reference points 
and luminosities are used. 

parameters are and Amg^, where Arrig^ is assumed to be exactly known. For both MINOS 
and the CNGS experiments, we use the background uncertainties given in Refs. [12,44,45], 
i.e., 10% for MINOS and 5% for ICARUS and OPERA. Then we simulate data for each value 
of Am|^ with 6'i3 = 0, and fit these data with as the only free parameter. This simplified 
procedure leads to a limit similar to (sin^ 2^i3)eff, which represents the ability to identify a 
signal (but not to extract sin^26'i3), and is similar to a simple estimate of S/ yS + B. 

In Figure ^1 we compare the discussed effective sin^ 2^13 limit to the ones of the experimen- 
tal collaborations. The solid black curves represent our results, whereas the dashed gray 
curves are taken from Refs. [12,44,45]. Within the Super-Kamiokande allowed atmospheric 
region, our simulation is in very good agreement with the results of the different collabo- 
rations. The slight deviation in the OPERA curve at large values of Amg^ comes from the 
efficiencies in Ref. [45], since they are only given as energy-integrated quantities. Thus, it 
is not possible to fully reproduce the energy dependence of the events. This effect becomes 
stronger if the oscillation maximum is shifted to higher energies compared to the reference 
point used in Ref. [45]. However, the influence on our results is marginal, since OPERA does 
not contribute significantly to the ^13 sensitivity. 

Note that the numbers given in Table El and the results shown in Figure ^1 do not allow a 
comparison of the different experiment performances, because the reference points used for 
the calculation of Table El are rather different, and so are the integrated luminosities. For a 
comparison on equal footing we refer to Table (HI and the discussion in Section 

B Simulation of the reactor experiments 

For our simulation of the reactor neutrino experiments, we closely follow our previous 
work in Ref. [20]. For the analyses presented here, we assume a near detector baseline 
of Lnd = 0.15 km and we consider two options for the far detector baseline: L-pT> = 1.05 km 
corresponds to the baseline at the CHOOZ site, whereas a baseline LpD = 1.7 km is close to 
the options considered for several other sites. We always fix the number of reactor neutrino 
events in the near detector to 2.94 ■ 10^, which implies that it has the same size as the 
far detector at Lpo = 1.05 km with 60 000 events (assuming the same efficiencies in both 
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(sin^ 2^13 )eff sensitivity limit (sin^ 2^3 )eff sensitivity limit (sin^2^i3)eff sensitivity limit 



Figure 10: The effective sin^ 2^13 sensitivity limit (as discussed in the main text) at 90% CL as a 
function of the true value of Aml^ for the conventional beam experiments. The exposure is lOkty for 
MINOS, 20kty for ICARUS and a nominal running time of five years for OPERA, corresponding to an 
exposure of 8.25 kty. The dashed gray curves come from the collaborations of the individual experiments 
and are taken from Ref. [44] for MINOS, Ref. [12] for ICARUS, and Ref. [45] for OPERA. The black curves 
are obtained for sin^ 2^23 = 1 and Am^i = sin^ 2^12 = Sep = for systematics only, i.e., correlations and 
degeneracies are not included, as for the dashed curves. 



detectors^). We allow an uncertainty of the overall reactor neutrino flux normalization of 
Cabs = 2.5%. For the normalization error between the two detectors, we use a typical value 
of (7rei = 0.6%. As showH in Ref. [20], this roughly corresponds to an effective normalization 
error of iTnorm — 0.8%. The total range for the visible energy Eyis — Ep — A + rUe (where 
A is the neutron-proton mass difference, and rUe is the electron mass) from 0.5 MeV to 
9.2 MeV is divided into 31 bins. Furthermore, we assume a Gaussian energy resolution with 
(Tres = 5%/ \/ £^vis [MeV] . We remark that our results do not change if a smaller bin width 
is chosen, as it would be allowed by the good energy resolution and the large number of 
events. Furthermore, we take into account an uncertainty on the energy scale calibration 
(7cai = 0.5%, and an uncertainty on the expected energy spectrum shape Ushapc = 2%. which 
we assume to be uncorrelated between the energy bins, but fully correlated between the 
corresponding bins in near and far detectors (see Ref. [20] for details). 

In addition to the analysis as performed in Ref. [20], we have investigated in greater detail 
the impact of a background for a reactor experiment of the Double-Chooz type. We take 
into account four different background sources with known shape: 

• A background from spallation neutrons coming from muons in the rock close to the 
detector. This background can be assumed to be flat as a function of energy to a first 
approximation (see, e.g., Figure 48 of Ref. [10]). 

• A background from accidental events. A 7 from radioactivity is followed by a second 
random 7 with more than 6 MeV faking a neutron signal. Those events are important 
for low energies. 

^Due to a higher background rate in the near detector, there will more dead time than in the far detector. 
This reduces the number of events in the near detector roughly by a factor of two. We have checked that 
this has a very small impact on the final sensitivity. 
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• Two correlated backgrounds from cosmogenic ^Li and ^He nuclei. Both are created 
by through-going muons and give /5-spectra with end points of 13.6 and 10.6 MeV, 
respectively. 

In Table IHl we give for each background a realistic estimate [52] for the expected number 
of events in near and far detectors relative to the number of reactor neutrino events. Since 
the near detector will have less rock overburden than the far detector, more background 
events are expected. They are, however, compensated by the much larger number of reactor 
neutrino events in the near detector. Therefore, we assume to a first approximation the same 
relative sizes for the backgrounds in near and far detector. In Figure ITT| the spectral shape 
of these backgrounds is shown. We assume that these shapes are exactly known, however, 
the overall normalization of each of the four background components in the two detectors 
is allowed to fluctuate independently with an error of ctbg = 50%. In addition to these 
backgrounds with known shape, we include a background from an unidentified source with 
a bin-to-bin uncorrelated error of 50%. We assume a total number of background events of 
0.5% of the reactor signal in both detectors, and a fiat energy shape for this background. 

As a general trend, we find that backgrounds with known shape do not significantly affect 
the sensitivity. This holds independently of the integrated luminosity. Even increasing 
the numbers given in Table IHl by a factor five does not change the picture. This behavior 
can be understood in terms of Figure where we show the signal (the spectrum without 
oscillation minus the spectrum for sin^ 29i3 = 0.05) and its statistical error compared to 
the various background spectra. For illustration, background levels significantly larger than 
in Table IHl are assumed. From this figure, it is obvious that the spectral shape of the 
signal is very different from that of all of the background components. Already at modest 
luminosities, such as for Double-Chooz, enough spectral information is available to determine 
the backgrounds with sufficient accuracy, which means that it is not possible to fake the 



Background type Spectral shape 


BG/Reactor events 


0"BG 


Backgrounds with known shape 


Spallation neutrons Flat 


0.4% 


50% 


Accidentals Low energies 


0.2% 


50% 


Cosmogenic ^Li /9-spectrum (end point 13.6 MeV) 


0.2% 


50% 


Cosmogenic ^He /3-spectrum (end point 10.6 MeV) 


0.2% 


50% 


Bin-to-bin correlated BG total: 


1.0% 




Bin-to-bin uncorrelated background 


Unknown source Flat 


0.5% 


50% 



Table 6: Backgrounds included in our reactor experiment analysis. For each background source, the 
column "BG/Reactor events" refers to the total number of background events in the energy range between 
0.5 and 9.2 MeV relative to the total number of reactor neutrino events for no oscillations [52]. We assume 
the same magnitudes of the backgrounds relative to the total events in the near and far detectors. For the 
backgrounds with known shape, ctbg is the uncertainty of the overall normalization. For the uncorrelated 
background, ctbg is the error on the number of events in each bin, which is uncorrelated between different 
bins (31 bins). All backgrounds are uncorrelated between the two detectors. 
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Figure 11: Energy spectrum of the backgrounds from spallation neutrons, accidentals, and cosmogenic 
^Li and '^He. The percentage given for each curve corresponds to the total number of background events 
relative to the total number of reactor neutrino events for no oscillations. Also shown is the total background 
spectrum. The curve labeled "signal" corresponds to A''(i?vis; sin^ 2^13 = 0) — iV(i?vis; sin^ 2^13 = 0.05), 
where N{Evis; sin^ 2^13) is the energy spectrum for given sin^ 2^13. The shaded region is the statistical error 
band at la, i.e., ±^7V(£'vis; sin^ 2^13 = 0). Note that the absolute normalizations of the backgrounds are 
exaggerated in this figure. 



signal within the statistical error by fluctuations of the background components. In contrast, 
we flnd that the bin-to-bin uncorrelated background only plays a minor role for experiments 
of the size of Double-Chooz, whereas it becomes important for large experiments, such as 
Reactor-ll. In the latter case, a background level of 0.5% with a bin-to-bin uncorrelated error 
larger than about 20% would signiflcantly affect the sensitivity. We conclude that for large 
reactor experiments, the shape of the expected background has to be well under control.'' 

In Table |71 we summarize the experimental parameters which we use for the two setups 
Double-Chooz and Reactor-ll in the main text of this study. For the Double-Chooz setup, we 
stick closely to the configuration discussed in Ref. [23]. We have checked that the effect of 
the slightly different distances (1.0 km and 1.1 km) of the far detector location from the two 
different reactor cores at the CHOOZ site is very small. Hence it is a good approximation 
to consider the average baseline of 1.05 km for both cores. In order to obtain robust results, 
we include all systematical errors as well as backgrounds. 

The large reactor experiment Reactor-ll corresponds to the same setup as already used in 
Ref. [20]. It represents an ideal configuration without backgrounds and any systematical 
errors beyond overall normalization errors. This setup has been chosen to illustrate how 

^Note that the above statements quantitatively depend to some extent on the chosen number of bins. 
For example, a bin-to-bin uncorrelated error of 20% for 10 bins has, in general, a different impact than such 
an error for 60 bins. This can be understood by the fact that for a large number of bins, the simultaneous 
fluctuation of neighboring bins becomes unlikely. The same argument holds for the flux shape uncertainty 

<^shape- 
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Double-Chooz 


Reactor-1 1 


TjiirninoRitv C 


288t GWv 


8000 1 GWv 


IMiimVjpr of pvpnts in IMTO 


2.94 • 10^ 


2.94 ■ lO*" 


IMiimliPT' of pvpnts in KTO 


6 ■ 10^ 


6.36 ■ 10^ 


Near detector baseline 


0.15km 


0.15 km 


Far detector baseline 


1.05 km 


1.70 km 


Fnprp'v rPKolntion 




5%/V^[MeV] 


Visiblp pnprEfv rantrp 


0.5 - 9.2 MeV (31 bins) 


Individual detector normalization 


(T,el = 0.6% 


drcl = 0.6% 


Flux normalization 


O-abs = 2.5% 


O-abs = 2.5% 


Flux shape uncertainty 


Cshape = 2% 


''^shape 


Energy scale error 


fXcal = 0.5% 


0"cal = 


Backgrounds included 


Yes 


No 



Table 7: Characteristics of the two reactor experiments Double-Chooz and Reactor-ll. 

an optimal reactor experiment would fit into the general picture of the next ten years of 
oscillation physics and, to obtain CP-complementary information. Several proposals which 
are close to our Reactor-ll setup are currently discussed [16]. 

C The definition of the sin^ 2^13 sensitivity hmit 

In this appendix, we discuss the definition of the sin^ 2^13 sensitivity limit. Although this 
definition is very general, we mainly focus on the v^. appearance channel, since one 

has to deal extensively with parameter correlations and degenerate solutions in this case. 
Let us first define our sin^ 26'i3 sensitivity and then discuss its properties. 

Definition 1 We define the sin^ 20i3 sensitivity limit as the largest value o/sin^26'i3, 
which fits the true value sin^ 2^13 = at the chosen confidence level. The largest value of 
sin^ 26'i3 is obtained from the projections of all ( disconnected) fit manifolds (best- fit manifold 
and degeneracies) onto the siv? 26i^-axis. 

Since for future experiments no data are available, one has to simulate data by calculating a 
"reference rate vector" for a fixed set of "true" parameter values. In general, the experiment 
performance depends on the chosen set of true parameter values, and it is interesting to 
discuss this dependency in many cases. It is especially relevant for the true values of ^m\^ 
and Amg^, which we usually choose within their currently allowed ranges. According to 
Definition we choose the true value sin^ 26'i3 = to calculate the sin^ 26'i3 limit, since 
we are interested in the bound on sin^ 2^13 if no positive signal is observed. Moreover, this 
choice has the following advantages: 

• Since for sin^ 2^13 = the phase ^cp becomes unphysical, the sensitivity limit will be 

independent of the true value of6cp- 



31 



• For sin^ 2^13 = 0, the reference rate vectors for the normal and the inverted mass 
hierarchies are approximately equal, which implies that the sensitivity limit hardly 
depends on the true sign of Aml^ (see also the discussion related to Figure [T21 later). 

Once the reference rate vector has been obtained, the fit manifold in the six-dimensional 
space of the oscillation parameters is given by the requirement Ax^ < CL {e.g., at the 
90% confidence level, we have CL = 2.71 for 1 d.o.f.). In addition to the allowed region 
which contains the best-fit point ("best-fit manifold"), one or more disconnected regions 
("degenerate solutions") will exist, and each of them may have a rather complicated shape 
in the six-dimensional space ("correlations"). The final sensitivity is given by the largest 
value of sin^ 26'i3 which fits sin^ 2^13 = 0. It is obtained by projecting all these disconnected 
fit regions onto the sin^ 26'13-axis, where the projection takes into account the correlations. 
Hence, this procedure provides a straightforward method to take into account correlations 
and degeneracies. Thus, for the case of the z/^ — > Ug appearance channel, our definition of the 
sin^ 2^13 sensitivity limit includes the intrinsic structure of Eq. (^J. This equation reflects 
that an appearance experiment is only sensitive to a particular combination of parameters. 
The projection onto the sin^ 2^13-axis takes into account that all the other parameters can be 
only measured with a certain accuracy by the experiment itself. Moreover, in complicated 
cases {e.g., for neutrino factories) local minima may appear in the projection of the x^- 
function onto the sin^ 2^13-axis, and the x^-function can intersect the chosen confidence level 
multiple times. In this case, we choose by definition the rightmost of these intersections. 
Hence, the sensitivity limit, as defined above, refers to the potential of an experiment (or 
combination of experiments) to extract the value of the parameter sin^ 26'i3 from Eq. (|T|) 
convolved with all the simulation information. 

In this study, we compare the sin^ 2^13 sensitivity limit of a given experiment to a so- 
called (sin^ 26*13)0^ sensitivity limit in some cases. The (sin^ 26'i3)eff sensitivity hmit roughly 
corresponds to the potential of a given experiment to observe a positive signal, which is 
parameterized by some (unphysical) mixing parameter (sin^ 2^i3)efr: 

Definition 2 The (sin^ 20i3)eff sensitivity limit is defined as the sensitivity limit from 
statistics and systematics only which is computed for 6cp = by fixing all other oscillation 
parameters to their true values. 

In order to illustrate the impact of systematics, correlations, and degeneracies, we often 
use "bar charts" (see, for example. Figures El and E)) , where the final sin^2^i3 sensitivity is 
obtained by successively switching on systematics, correlations, and degeneracies. In these 
bar charts, the statistics-only sin^ 2^13 sensitivity (left edge of the bar) is computed for all 
oscillation parameters fixed and Sqp = 0, the statistics+systematics sensitivity limit cor- 
responds to the (sin^ 2^i3)efr sensitivity limit, the statistics+systematics+correlations limit 
corresponds to the sensitivity limit for the best-fit manifold only (no degenerate solutions 
included) , and the final sensitivity limit (right edge of the bar) corresponds to Definition ^ 

In the following, we illustrate in greater detail how the sin^ 2^13 limit is obtained, how the 
bar charts are constructed, and how the sin^ 26*13 and (sin^ 2^i3)efr limits are related to each 
other at the example of the JPARC-SK experiment. We focus mainly on the sgn(Am3^)- 
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Inverted hierarchy 




0.01 0.02 0.03 

sin^29i3 sensitivity limit 



Figure 12: The as function of sin^ 26*13 for JPARC-SK. For the true values of the oscillation 
parameters, we choose the current best-fit values from Eq. ©, sin^ 2^13 = 0, ^cp = (for curves without 
correlations only), and normal (upper plot) or inverted (lower plot) mass hierarchies. The solid curves in 
each plot are obtained by fitting with the same mass hierarchy as has been used to calculate the reference 
rate vector ( "right-sign" ) , whereas for the dashed curves the wrong mass hierarchy has been used ( "wrong- 
sign"). Within each group of solid or dashed curves, the left curve determines the statistics-only limit, the 
middle curve the statistics-|-systematics limit, and the right curve the statistics-|-systematics-|-correlations 
limit. Note that the wrong-sign minimum has not exactly the same position in parameter space as the 
original minimum. 

degeneracy and the correlation between ^13 and 5cP) wliicli is of particular relevance for the 
appearance channel at superbeams. 

In Figure ^1 the is shown as a function of sin^ 26'i3 for the "right-sign" and "wrong- 
sign" solutions, where in the upper (lower) plot the normal (inverted) hierarchy has been 
chosen to calculate the reference rate vector. The right-sign solution is obtained by fitting 
with the same sign of Amgj^ as the reference rate vector has been calculated with, z.e., 
the "right" neutrino mass hierarchy is used, whereas the wrong-sign solution is obtained 
by fitting with the opposite sign of Amg^^, z.e., the "wrong" mass hierarchy is used. The 
different curves in each group with the same curve style correspond, from the left to the 
right, to the statistics-only, statistics-|-systematics, and statistics-|-systematics-|-correlations 
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Figure 13: The 90% CL fit manifold (1 d.oi.) in the sin^ 26'i3-(5cp-plane for JPARC-SK. For the true 
values of the oscillation parameters, we choose the current best-fit values from Eq. ((2Jl and sin^ 26'i3 = 0. The 
different curves correspond to various sections (un-displayed oscillation parameters fixed) and projections 
(minimized over un-displayed oscillation parameters) as described in the plot legend. The bars demonstrate 
the individual contributions to the final sin^ 29i3 sensitivity limit. Note that for the sgn(Am|]^)-degenerate 
solution, we only show the final projection. 

sensitivity limits, where these hmits are obtained from the intersection of the with the 
Ax^ = 2.71 hne. The bar charts are constructed from the corresponding curves, as one can 
easily read off the figure. 

Comparing the normal and inverted mass hierarchy plots in Figure ^1 one can observe a 
symmetry between the right- and wrong-sign solutions: The curves for the normal mass 
hierarchy and Am|]^ > (right sign) are very similar to the ones of the inverted mass 
hierarchy and Am|]^ > (wrong sign). This can be understood in terms of the identical 
appearance rate vectors for the normal and inverted mass hierarchies for the true value of 
sin^ 2^13 = 0. However, since the role of the Amg^^ > curves is different for the normal 
and inverted mass hierarchies, i.e., they either correspond to the best-fit manifold (right 
sign) or the sgn(Am|^)-degeneracy (wrong sign), the bar charts are, by definition, very 
different, since they are originally determined by the best-fit solution. However, one can 
easily see that the final sensitivity limit does not depend on the mass hierarchy [26]. This 
property comes from the fact that the degeneracy part does not contribute to the final 
sensitivity if the best-fit sin^ 2^13 sensitivity is already worse than the degenerate solution 
sensitivity. Since there is hardly a difference between final sensitivity limits for the different 
mass hierarchies, we usually show the normal mass hierarchy sensitivity limit. In fact, there 
is a small difference between the final sensitivity limits for the different mass hierarchies, 
which mainly comes from the disappearance channels. 

Let us now illustrate the impact of the correlation between sin^ 26*13 and 6cp. Therefore, 
we show in Figure ^1 the fit manifold in the sin^ 2^i3-5cp-plane. The sin^ 2^13 sensitiv- 



34 



ity limit is again obtained from the projection onto the sin^ 2^13-axis. In Figure ^1 the 
individual contributions to the bar chart are illustrated by showing different sections (un- 
displayed oscillation parameters fixed, i.e., no correlations) and projections (minimized over 
un-displayed oscillation parameters, i.e., they include correlations) of the fit manifold. One 
can see that both edges of the leftmost (blue) bar are computed for 6cp = 0. This illustrates 
that if 6cp and all the other oscillation parameters except sin^ 2^13 are fixed at the true val- 
ues, much stronger bounds on sin^ 26*13 can be obtained, corresponding to the (sin^ 26*13)65 
limit. The limit gets considerably weaker if the x^-function is minimized over Sqp, as well 
as all oscillation parameters which are not shown, which leads to the "correlation bar" . In 
fact, from Figure IT^ one can see that the largest part of the correlation bar comes from the 
correlation with 6cp [53], and only the small difference between the dark dashed and the 
light solid curves comes from the correlation with the other oscillation parameters. The final 
sensitivity limit is then obtained as the maximum value of sin 2^13 which fits sin^ 2^13 = 
including all degenerate solutions. 
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